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We discuss solar system constraints on f{Q) gravity models, where / is a function of the Gauss- 
Bonnet term Q. We focus on cosmologically viable f{Q) models that can be responsible for late-time 
cosmic acceleration. These models generally give rise to corrections of the form e(r/r s ) p to the 
' vacuum Schwarzschild solution, where e = H^rl <C 1, r s is the Schwarzschild radius of Sun, and 

, H, is the Hubble parameter today. We generally estimate the strength of modifications to General 

■ Relativity in order to confront models with a number of experiments such as the deflection of light 

' and the perihelion shift. We show that cosmologically viable f{G) models can be consistent with 

solar system constraints for a wide range of model parameters. 

I. INTRODUCTION 

The modified gravity paradigm has been extensively studied over the last few years as a way to make gravity 
responsible for the observed acceleration of the universe at large scales [l[ . These models are attractive in the sense 
that cosmic acceleration can be realized without recourse to a dark energy matter component. Unlike the cosmological 
O i' constant scenario, they generally give rise to a dynamical equation of state of dark energy that varies in time. 

Among these modifications of gravity, the so called f(R) theory has become popular and started to be a new 
branch of cosmology Q. In this theory the Lagrangian density / is function of the Ricci scalar R. The f(R) theory 
in the metric formalism can be proven to be classically equivalent to a class of Brans-Dickc theory with a Brans- 
Dicke parameter wbd = [jj. A scalar-field degree of freedom, called scalaron 0], can freely propagate to mediate a 
^ ' fifth force, unless the scalaron mass is heavy in the region of high density. For the compatibility with local gravity 
experiments the f(R) theory needs to approach the Lagrangian density f(R) — R — 2A in General Relativity (GR) 
for large values of R much larger than the present cosmological Ricci scalar i?o @ ■ Meanwhile non-trivial deviation 
from the ACDM model can arise for R of the order of R j& 0, 0, @, G3 [111 ]. This leads to a number of interesting 
observational signatures such as the modifications to the matter power spectrum [l2| as well as to the weak lensing 
spectrum [l3j . 

It is also possible to take into account a Gauss-Bonnet (GB) scalar Q that is a combination of the Ricci scalar R, 
the Ricci tensor R^ v , and the Riemann tensor Rp, va j3 [Til. [ill]. This GB scalar, together with R, belongs to an infinite 
class of curvature invariants, the Lovelock scalars, which have the property that they do not introduce derivatives 
terms higher than two into the equations of motions for the metric tensor. Among these scalars, R and Q are the 
only ones that do not identically vanish in four dimensions (apart from the cosmological constant). However the term 
\J—g G is a total derivative in four dimensions, where g is the determinant of the metric tensor. This means that the 
only way for the GB term to contribute to the equations of motion is to couple it to something else, e.g., a scalar 
field 4> with the coupling of the form F(<j>)Q fl6| . This kind of coupling is present in the low energy effective action of 
string-theory (l7j . due to the presence of dilaton-graviton mixing terms. 

The dynamics of dark energy based on the dilatonic coupling F(<p) oc with the exponential potential V(cf>) cx e~ x< ^ 
has been studied by a number of authors [HI, 0, HH H3] . While it is possible to realize a scaling matter era 0, [2Cjj 
followed by a late-time cosmic acceleration, the region of viable parameter space to satisfy several observational 
constraints is restricted to be very small Um. It was also shown in Refs. [2(| [23| that tensor perturbations tend to 
exhibit negative instabilities if the GB term is responsible for cosmic acceleration. Moreover, in such models, the 
energy fraction of the GB term needs to be strongly suppressed for the compatibility with local gravity experiments 
[24[, which is at odds with the requirement of cosmic acceleration induced by the GB term. 

There is another class of modified gravity models in which the Lagrangian density is described by R+f(G) (so-called 
"f(G) gravity"), where / is function in terms of the GB term Q j25[ . Unlike f(R) gravity, this theory does not have 
an action in the Einstein frame with a standard kinetic term of a scalar-field degree of freedom. The conditions for 
cosmological viabilities of f(G) gravity have been studied in Refs. [26|, [2?], [28| (see also Refs. j29l. [30|). Li et al. [26j 
showed that the condition < H 6 f gg <C 1 (where f t gg = d 2 f/dG 2 ) needs to be fulfilled in order to keep cosmological 
perturbations under control. In Rcf. [271 ] the same condition has been derived to ensure the stability of a late-time de 
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Sitter solution as well as the existence of standard radiation and matter dominated epochs. In particular the stability 
of the de Sitter point requires the condition < Hffgg < 1/384, where Hi is the Hubble parameter at this point. In 
order to remove possible singularities in the cosmic expansion history the second derivative fgg should not change 
sign, i.e. f t gg > for all G, together with the condition that f^gg — > as \G\ — ► oo. This removes the presence of 
unstable modes during the whole past evolution of the universe. A number of cosmologically viable models satisfying 
these requirements have been proposed in Ref. [27j ]. 

In this paper we will study the property of f(Q) gravity on solar system scales and investigate whether cosmologically 
viable f(G) models can be consistent with solar system constraints. We first find approximate vacuum solutions 
for these models through an iterative method introduced in Ref. [3l| . The solutions look like corrections to the 
Schwarzschild solution, where the corrections are typically in the form of positive powers in the ratio r/r si where r is 
the distance from the center of compact object and r s is the Schwarzschild radius. This behavior of positive powers 
is similar to the typical correction that the cosmological constant gives to the Schwarzschild solution (oc (r/r s ) 2 ). In 
the case of the Schwarzschild-de Sitter solution the metric in the solar system is dominated by the term r s /r, but one 
can put upper bounds on the value of the cosmological constant such that its contribution is allowed by experimental 
data. We follow a similar procedure in order to constrain the values of model parameters in f{G) gravity. We will 
show that cosmologically viable f(G) models can satisfy solar system constraints for a wide range of parameter space. 

In section II wc briefly review cosmologically viable f(G) models. In section III we describe the method used to 
find approximate spherically symmetric solutions of the Einstein equations. In section IV we discuss solar system 
constraints in the presence of positive powers of the radius correction to the Schwarzschild solution. In sections V, 
VI, and VII we apply the constraints to a number of f(G) models. In section VIII we report our conclusions. 

II. COSMOLOGICALLY VIABLE f(Q) MODELS 

Let us first briefly review cosmologically viable f{G) models proposed in Ref. [27j • The action describing this theory 
is given by 



$m\9nvi ^m) j 



(1) 



where R is a Ricci scalar, Q = R 2 — AR^ V R^ V + R^apR 11 "" 13 is a Gauss-Bonnet (GB) term, and S m is a matter action 
that depends on a spacetime metric g^ v and matter fields ^> m - We use the unit M p i = 1/ \Z8itGn = 1, but we restore 
the reduced Planck mass M p \ and the gravitational constant Gjy if required. 

The late-time cosmic acceleration can be realized by the presence of a de Sitter (dS) point satisfying the condition 
Gif,g(Gi) — f(Gi), where Hi and Gi are the Hubble parameter and the GB term at the dS point respectively. 
The stability of the dS point demands the condition < Hflgg(Hi) < 1/384 |27|. The GB term, Q = 24H 2 (H 2 +H), 
changes sign from negative to positive during the transition from the matter era to the accelerated epoch. For the 
existence of standard radiation and matter eras we require that f_gg > for G < Gi and that f t gg approaches +0 in 
the limit \G\ ~ * oo. We also need the regularities of f(Q) and its derivatives f t g, f t gg- The following two models can 
satisfy these conditions (27| : 



3H 2 



(A) /(^)=A-^=arctanf|- 

(B) f(G) = \-j= arctan 
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aXy G 



(2) 
(3) 



where a, A and G* are positive constants. Note that G* roughly corresponds to the scale Hf for a and A of the 
order of unity, where H* is the Hubble parameter today. The second derivative of / with respect to G is f,gg = 

X/[G* /2 (l + G 2 IGl)\ and f.gg = 2A/[g* 3/2 (l + G 2 /G 2 ) 2 } for the models (A) and (B) respectively, so that f t gg > for 
A > 0. 

In the region of high density where local gravity experiments are carried out (\G\ S> G*), the above models have the 
following asymptotic behavior 



(A) f(G) 

(B) f(G) 



1 i G 

2 vu: 
2 vg: 



i— 1 i— (G 2 
(a+l)Av^- 2 X ^ hx [jj2 



(a + l)AV^ + 



WG^ql 

3 G 2 



XVGZGl 
6 G 2 



(4) 
(5) 



3 



The first terms in Eqs. (|4]) and §5§ arc linear in Q so that they do not give rise to any contribution to the Einstein 
equation, whereas the second terms contribute to the field equation as a cosmological constant. The other terms in 
Eqs. (j4]) and correspond to the corrections to the ACDM model. The difference between the models (A) and (B) 
is that the former has a logarithmic correction that mildly increases with the growth of \Q\. Note that the viable f(R) 
models such as (i) f(R) = R : - XR*(R/R*) 2n /{(R/R t ) 2n + 1] and (ii) f(R) = R - AiJ»[l - (1 + R 2 /Rl)- n ] (n > 0) 
behave as f(R) — R — XR* + \R*(R*/R) 2n in the region of high density (R 3> R*)- This asymptotic form is similar 
to the model (B) given above by replacing R for Q. 

In the following we shall study solar-system constraints on cosmologically viable f(G) models. Before doing so, it is 
worth mentioning the difference between f(R) and f(Q) theories. If we consider a spherically symmetric background, 
the Schwarzschild vacuum solution corresponds the vanishing Ricci scalar (R = 0). In the presence of non-relativistic 
matter, R approximately equals to the matter density /9 m /M 2 j for viable f(R) models Q EH- Then the term 
{R st /R) 2n is roughly of the order of (p c /p m ) 2n , where p c is the cosmological density today. The ratio (p c / p m ) 2n 
becomes much smaller than 1 for n > in the region of high density (p m 3> p c ) so that one has (R*/R) 2n <C 1. In 
the presence of non-relativistic matter the chameleon mechanism [32| can be at work for the f(R) models that have 
the asymptotic form f(R) — R — XR* + XR*(R*/R) 2n in the region R which allows the possibility for the 

consistency with local gravity tests. In fact it was shown in Ref. [33[ that these models can satisfy solar system and 
equivalence principle constraints for n > 0.9. 

On the contrary one has Q = R llva pR' lva0 = 12 r 2 /r 6 on the vacuum Schwarzschild solution, where r s = 2GnMq 
is the Schwarzschild radius and M is the mass of the star. Since Q does not vanish even in the vacuum, the term 
such as (G 2 /G 2 ) n (n > 0) can be much smaller than 1 even in the absence of non-relativistic matter. If matter is 
present, this will give rise to the contribution of the order of R 2 w (8ttG n p m ) 2 to the GB term. The ratio of the 
matter contribution to the vacuum GB value = 12 r 2 /r 6 = 48(GatM ) 2 /r 6 can be estimated as 



R 2 ^ {^) 2 p 2 m r* 



g(°) 48 M 2 



(6) 



As long as s <C 1, we can neglect the matter contribution to the GB term. 

At the surface of Sun (radius r Q = 6.96 x 10 10 cm= 3.53 x 10 24 GeV _1 and mass M = 1.99 x 10 33 g = 1.12 x 
10 57 GeV), the density p m drops down rapidly from the order p m w 10 -2 g/cm 3 to the order p m w 10~ 16 g/cm 3 . If we 
take the value p m = 10 _2 g/cm 3 we have s f» 4 x 10 -5 (where we have used lg/cm 3 = 4.31 x 10~ 18 GeV 4 ). Taking 
the value p m = 10~ 16 g/cm 3 leads to a much smaller ratio: s«4x 10 -33 . The matter density approaches a constant 
value p m ~ 10 - g/cm 3 around the distance r = 10 3 tq from the center of Sun. Even at this distance we have 
s«4x 10 -31 , which means that the matter contribution to the GB term can be completely neglected in the solar 
system we are interested in. If we use the value p m « 10 -24 g/cm 3 , s exceeds the order of 1 for the distance r > 10 8 r Q . 
However this is out of the region where solar system experiments are concerned. Moreover the Schwarzschild solution 
has no meaning far away from the star where other contributions can arise, i.e. other close stars, the mean field of 
the galaxy, and so on. From the above discussion we find that vacuum solutions can be used when we discuss solar 
system constraints on f(Q) gravity. 



III. EXPANSION AROUND THE SPHERICALLY SYMMETRIC SPACETIME 

The cosmologically viable f(G) models ^ and ([3]) will consist of a numerical factor of order \fQl times a dimen- 
sionlcss function (because f(G) has the dimension of [mass] 2 ). Meanwhile the typical value of the GB term on the 
vacuum Schwarzschild solution is given by 

g s = U/rt. (7) 
When we discuss solar system constraints, it is convenient to define the following dimensionless ratio 




Since s/G* is of the order of the squared of the present Hubble parameter H* , the parameter e is approximately given 
by e ~ (i/*?* s ) 2 /(2\/3). In the subsequent sections we shall discuss the case of Sun with the Schwarzschild radius 
r s = 2.95 x 10 3 m. Using the value w 70 km sec -1 Mpc -1 , the parameter e for Sun is approximately given by 



e w 10~ 46 



(9) 
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The model (|2|) can be written in the form f(Q) — e \\\/Glg{x) — aXy/G7\ , where g(x) = xarctanx — (1/2) ln(l + x 2 ) 
and x = G/G* = G/(G S £ 2 )- Hence the function / can be replaced by the form / = ef. 
The equations of motion in the vacuum can be written as 

G% + e = , (10) 

where G^ y is the Einstein tensor, and 

= 8 [R^a + R pv gan - Rpagvv - R^gap + Rpo-g^p + R{g^g<j P - svsw>)/2] VV7, g + (Gf,g - f)gp V - (11) 

In general these theories will have vacuum solutions, which we wish to study here. Although exact solutions are not 
always found analytically, it is possible to obtain approximate solutions which reproduce the real ones very well at 
least on some scales. In fact, since e <C 1, we can try to use the iterative method introduced in Ref. 
We look for static spherical symmetric solutions of the kind 

dr 2 

ds 2 = -A(r,e)dt 2 + — + r 2 (d(9 2 + sin 2 dip 2 ) , (12) 

B{r,e) 

where the functions A and B are Taylor expanded in the form: 

A(r, e) = A (r) + ^i(r)e + A 2 (r) e 2 + . . . , (13) 
B(r,e) = B (r)+B 1 (r)e + B 2 (r)s 2 + ... . (14) 

Using this expansion for A and B, both G M „(r, e) and £ M „(r, e) can be expanded order by order in e. At lowest level, 
the equations of motion correspond to 

G» v ^=0, (15) 

which leads to the usual Schwarzschild solution Bq = Aq = 1 — r s /r. At first order one has 

e[G» v ( 1) {A 1 ,B 1 ,A ,B )+V v ( -°\A ,B )} = 0, (16) 

which reduces to ordinary differential equations for Ai,B±, where Aq, Bq are the Schwarzschild terms found previously. 
This method can be reiterated to get the coefficients A2,B2, as well as all the other higher-order terms. 

It should be noted that for general approach to the modification of gravity, including quintessence, there is no more 
analogue of the Birkhoff theorem regarding the unicity of the solution for a spherically symmetric vacuum solution. 
We can only impose to have a static metric. The bottom line is that the forms of A and B can be found at any 
order-e approximation, by solving the Einstein equations iteratively with respect to the small parameter 0<e<1. 

Suppose that we find such an iterative solution and write both A and B as power expansions of e. Then the iterative 
solution should have in general a radius of convergence, inside which each correction is larger than the next order e 
term. In this case we expect that the dominant term corresponds to the Schwarzschild contribution. Therefore in the 
small e limit we will have 

A = 1 - l/p + e ci (f + e 2 c 3 p m + 0(e 3 ) , (17) 
B = l-l/p + £c 2 pQ +e 2 c iP n +0(e 3 ), (18) 

where p = r/r s . In Sees. IVII and [VIII we will show that the correction terms in Eqs. (fTT)) and ([THfl in fact arises for 
the cosmologically viable f(G) models given in ^ and We will restrict ourselves to the region p ^> 1, which is 
generically satisfied in the solar system outside Sun. In other words this corresponds to the weak limit of the theory. 
We have introduced ci^.3,4 as constants whose values need to be bounded experimentally, and also p,q,m,n as the 
powers of p. In the following we will restrict our attention to the case p = q > and m = n > 0. Later we will see 
that this assumption is quite general for cosmologically viable f(G) models. 

In order to have a meaningful e expansion, one needs to verify that there exists a set of convergence, e.g., 

D = {p|p»l, sff^p- 1 }. (19) 



This implies that the expansions (|17[) and (|18p can be trusted only in this region. Outside it, one must find the 
solutions of Einstein equations, both analytically or numerically, without using the e expansion. This also implies 
that it does not make sense to study this solution in the limit p — > 00. However, even if the full solution is known, it 
does not have physical meaning for values of r where the solar system begins to feel other stars' contribution. This 
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situation is analogous to the Schwarzschild-dc Sitter solution, where the same solution cannot be trusted any more at 
distances a few parsecs away from Sun. In this region the solar system cannot be treated as an entity isolated from 
the rest of the galaxy. Of course, for larger p, the set D becomes smaller. However, this is not enough: one should 
also verify that in the same set the e 2 -term is much smaller than the e-one. This implies that in D one also requires 
to have 

D 2 = {p | p » 1, e P m « pO}CD. (20) 

If < m < 2p + 1 then this condition is verified. We will call a good e expansion for the one where this last condition 
is valid. In general we may have more complicated cases where there exists an order of expansion d such that Di C Dd 
for all i > 1 (and d might not be one), such that the expansion makes sense in Dd- If there is no such set, then the 
expansion does not make sense. If there is such a set, we can define D = Dd- 

If D 2 C D, then, in D (or, more in general, in Dd), it is safe to approximate the perturbative solution as 

A= l-l/p + eap" , B = 1 - l/p + ec 2 p p . (21) 

In the next section we shall study a number of solar system bounds for the metric (|21[) . 



IV. SOLAR SYSTEM CONSTRAINTS 



There are a number of solar system constraints on the deviation from General Relativity-such as (a) deflection 
of light, (b) Cassini experiment, (c) Perihelion shift, (d) retardation of light, and (e) gravitational redshift. In the 
following we discuss those constraints for a general metric in the form of (|2ip . 



A. Deflection of light 

The first constraint we discuss is the deflection of light. The Lagrangian for a photon moving in the 9 = ir/2 plane 
in the gravitational field of the metric (TT2|) is given by 

L = lAP- i-ETV 2 - §rV -0, (22) 

where a dot represents an affine parameter along the geodesies. There are two constants of motion, namely 

E = At and L = r 2 ip . (23) 

Then we find 

, 9 f E 2 B B\ 

The minimal distance ro can be defined such that r(ro) = 0, giving 

T 2 r? 

r%E 2 A(po) ' { 0) 

where p = r /r s . 

Integrating dip/dr = ip/r by using Eqs. ([23|) -([25 ]) . we obtain 




where p is used to avoid the confusion with the upper limit p of the integral. The Schwarzschild-de Sitter solution 
corresponds to c\ = c 2 and p = 2, in which case the integral (|2l?)) reduces to the standard GR contribution. Therefore 
the cosmological constant does not give any modification to GR for light bending [Hj . In standard GR, the integral is 
taken in the limit p — > oo. This is a mathematical extrapolation, as the metric does not hold in the whole spacetimc. 
Of course this property holds in our approach where the approximate metric is sensible only in D, the domain of 
convergence. Up to first order in e and in the domain £>, it is possible to approximate the integral as 

<p = ± [ P dp^J _ f P0 3f ^ ;+e f P dpF(p ,p), (27) 

J pa P V P (PO ~ 1) - PoKP ~ !) Jpo 
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where 

1 3/2 (f{c lP \ PQ - 1) + c 2 [p 3 (l - p ) + ( P - l)pl] }- Cl (p- l)pV +1 



F( Po ,p) 



PPo 



p 3 ( Po -i)- ppI + pI 



2p(p-l) 



(28) 



We can further simplify this expression by considering the limits p, po 3> 1: 



•/*> P\/P rZr pl Jp» 2p 2 (p + p )^p 2 -p 2 2 7 P0 p(p 2 -p2)3/2 

Taking the positive sign in front of each integral, the deviation angle in the region D is given by 

<&(p)=2 dp = + / dp—— + ep / dp — ^ 7T, (30) 

•/po P ^P " P§ -Vo P 2 (P + PoVp "Po ^Po P(P 2 -Po) 3/2 

in the limit that p 3> po- The first two contributions in Eq. (|29[) . corresponding to the GR ones, give 

#gr = — + 0(po/p) ■ (31) 
Po 

Meanwhile the e-contributions can be evaluated as 

— eci po + epo(ci - c 2 )ln(2p/p )+0(Po 1 ), (32) 
(ci-c 2 ) £ p p + 0(po 1 ), (33) 
|(c! - c 2 )ep p 2 + cep 3 [ln(2p/p ) - 1] + O(p x ), (34) 

^-^ e po p^ 1 + -At £ P? P P " 3 + G(pr 5 ) • (35) 
pi p o 

We have kept not only the dominant correction but also the next one. Of course the latter becomes important only 
for ci = c 2 , that is, when the first correction vanishes. Therefore we have evaluated the second (smaller) contribution 
only when the first vanishes, that is, when c\ = c 2 . In other words, we define c = c\ if c\ = c 2 . The experimental 
bound on <& e relative to f?GR is given by the Very Long Baseline Intcrferometry (VLBI), which combined measurements 
taken by different radio telescopes on Earth [35(. The experimental value i? cxp relative to the theoretical prediction 
has been constrained to be $ e xp/^theor = 1.0001 ± 0.0001. Therefore, in order that fl e does not affect the GR result, 
we impose 

< 10- 4 . (36) 



p = l, 




p = 2, 




P = 3, 


•&e 


P>4, 


# 6 



When p > 2 this condition translates into 



Recall that in the domain D we have 



2(p-l) 



ep^- 1 <KT 4 . (37) 



£PoP p - 1 «-^|«-, (38) 
Po P Po 

in which case the condition |37|) is satisfied for ci , c 2 , p of the order of unity. Hence the deflection of light always 
remains a small correction in the domain of convergence. 



B. Cassini experiment 

less et al. [36| showed that the contribution to the fractional frequency shift of a stable and coherent two-way radio 
signal (Earth-spacecraft-Earth) y, due to the metric of a gravitational theory (which possesses a weak field limit), is 
proportional to the deviation angle $ of light, and it is given by the formula 



(39) 
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where vq and v\ are transverse velocities of Earth and a spacecraft, and Iq and l\ are their distances from Sun. Since 
h 3> Iqi the £ contribution can be written as 

y E f» 2w E ai'th [<p(pO, PEexth) + <p(j>0, PCassini)] 

« «Earth °* _^ £ PO (Pcassi n i + ^Earth) » ( 40 ) 

where the approximate equality in the second line is valid for p > 2. Meanwhile, from Eq. pip , the GR contribution 
can be written as 

_ ., u Earth , . 1 

2/GR ~ 4 . (41) 

Po 

The e contribution needs to be a negligible correction to GR. In order to have y £ <C J/gr we require that 

epgpg^ini « 1 , ( 42 ) 

which is satisfied, as we have seen before, in the region D, of course if pcassini £ D. 

The signal due to GR detected by Cassini is y ~ 10~ 10 , within an experimental error of Ay oxp - 1CT 14 . Therefore 
the contribution y £ from the modifications of gravity needs to satisfy the condition y e < Ay cxp ~ 10 , or equivalently, 
ye/yGK < 10 -4 . If p > 2, this condition translates into 

l^^/Cini < IQ" 4 • (43) 

As we have shown in Eq. (J3SJ), this relation is satisfied in the domain of convergence. The constraint (|43p can be used 
to place experimental bounds on f(G) models later. 



C. Perihelion shift 



Let us proceed to constraints coming from the perihelion shift of the inner planets, whose GR value is in extremely 
good agreement with experimental data. We will follow a similar procedure to the one discussed in subsection II V Al 
The difference is that we now deal with the motion of a massive particle whose Lagrangian reads 

L=\Ai 2 -\B- 1 r 2 -\r 2 v 2 = \- (44) 
As in the previous case, the constants of motion are 

E — At and L = r 2 tp. (45) 
Defining u = r s /r, we can derive the following differential equation 

d 2 u rl r 2 E 2 B ( 1 dB 1 dA\ 1 dB , , . r 2 f dB 



' U -^ 2= ^a{b^-A^)-2^- {B - 1)u -^{ : ^ + 1] 
The flat-space solution can be obtained by setting the right hand side of Eq. (J3U) to be zero, as 



A 

2L 



a [l + *cos(0-0o)], (47) 



where < S < 1 is the eccentricity of the closed orbit. 
From Eq. (M it follows that 



The minimum distance po = fo/ r s satisfies 

r 2 E 2 ( 
A(p 



2L 2 v u/ \2p z 2L 2 



l -^ + "^\M + h)-k + ^- (49 > 
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Since both E and L are constants, Eq. (|49|) is an algebraic condition for po. Therefore, once we fix the orbit, that 
is v s , E and L, then we also fix r . The same solution will be valid at all times, therefore the orbit will always have 
the same condition for ro, i.e. the same perihelion and the same aphelion. In other words, the minimum (and the 
maximum as well) value for r will be unchanged at successive perihelia. Since r = also implies dr /dp = 0, the same 
initial conditions for the differential equation at the perihelia are identical at each perihelion and so the orbit repeats 
exactly, see Rindler's book [37l |. 

Expanding Eq. (|46[) at linear order in e, we find 



8 {c lP (P0 - 1) (i 2 r 7 2 + pi) + c 2 {L 2 r s - 2 [<p - 2)u 2 pg - p( Po 1)] + PP 2 }} ■ (50) 



— u 



dip 2 2L 2 2 2L 2 p% 

We evaluate the right hand side of Eq. ([50")) on the Newtonian solution u b . In this case, using Eq. (|4"T)) . we have 

L 2 1 



2 Po(l + S). (51) 



By doing so, Eq. (|50j) reduces to 

+ « - ^ = 2( niffZir Mi + 1)(po " + 2pa + 1} + C2(p " 2) * w v 

+ c 2 {p[S 2 (l - po) + 25 + 2p Q + 1] - 2p } + 2c 2 (p - 2)Sp cosy} . (52) 
In the limit po 3> 1, the right hand side can be simplified to 



d 2 u r 2 e[ PQ (l + 6)}P a 

oV + u - 2Z^ = ^TiW^TTj^ 12 ^ (1 + dH 

+ Po [cip (5 2 - 1) + ca(p - 2)5 2 cos 2 y - c 2 (p(<5 2 - 2) + 2) + 2c 2 (p - 2)<J cos y]} . (53) 



The second term on the right hand side of Eq. (|53j) is subdominant unless c\ — 0. In the following we will focus on 
the case c\ ^ 0. 

The solutions can be written down and studied for each p. As it happens in GR, the homogeneous solution can be 
described by the Newtonian solution with some periodic corrections, i.e. only dependent on cosy and higher harmonics 
(cos2y and so on). However there will be terms which have a secular impact on the orbit, and here we are looking 
exactly for such terms. The solution for Eq. (|53p , valid for any p, is given by 



r 2 3<5ysiny cip^P\{5 2 ) 

2L 2[l + 5 COS ^ ] + 2(1 + 5) 2 p 2 - S£piP+1) Vl-5 2 (l-S) P Sm ^ ai ' Ctan 



1-8 if 

TTs tan 2 



+ P 2 (6,<p), (54) 



where P\{x) is a polynomial such that Pi(0) = 1 with degree in x equal to the integer part of p/2, and Pi is periodic 
in y so that it does not have any secular contribution. The first contribution on the right hand side of Eq. (|54|) is 
the Newtonian solution, the second one is the standard GR secular correction, whereas the third one is the secular 
contribution due to the e-modification of the metric. Let us examine orbits with small eccentricity. Then at lowest 
order in 8 one finds that Eq. (f54"| reduces to 

r s r 3<5 y sin y 1 „ 

u ~ 7^2 l 1 + 6 cosy] + 2 (i + §)2 p 2 ~ -^6£c 1 p{p+ l)pgy siny, (55) 

where we have neglected a possible constant contribution from the arctan which, multiplied by sin y, reduces to a 
periodic term. 

The result (|55p coincides with another simpler method at lowest order in 5. One can expand Eq. (|53p in series of 
6, giving 

d 2 « r 2 . 1U „ 6 5 (cosy - l) + 3 „, r2 x 

— + u -^ = c l psp p [l-{p + l)8cos<p + p5} + !s ^5-^ +0{8 2 ). (56) 



The solution of this equation is 



38 y sin y l 



— ^2 [l + <5 cos y] H — -J £ £ ci p (p + l) p% (p simp + d\ cos(2y) + d 2 cos(3y) , (57) 

2Li 2po 2 
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where the coefficients are not needed to be known for our purpose. The second term of Eq. (|57|) coincides with 
the second term of Eq. (|55[) at lowest order in 5. The corrections are due to higher orders of 5, which are not included 
in Eq. (f57]). 

Using the relation (|51[) . the approximate solution (|55[) can be rewritten as 



where a 1) is defined by 



U ~ 2L? i 1 + 6 C0S (<P ~ v^)} , 



(58) 



(59) 



From Eq. (|58[) we find that in one orbit the angle between two perihelia is larger than 2ir approximately by 2na, or 



A<p - 2tt = 



3tt 



(1 + %o 



7T £ Cl p (jp +1)(1 + 5) p P 



P+l 



3 71 



(l + *)po 



l-i Cl ep(p+l)(l + 2<5K +2 



(60) 



The experimental bound on the shift Aip— 2tt for Earth, based on several thousands of optical observations of planetary 
positions, is [38[ 



Aip — 2tt — 5 ± 1 arcsec/century. 
For Mercury the bound is 43.1 ±0.1 arcsec/century. Since the GR contribution is given by 

one sets the modifications of gravity to contribute less than the experimental relative error, that is 



A Ve 



The correction term remains as such if 



l -\ Cl \ep{p+l)(l + 25) P l +2 <\. 



e P l +2 « 1 



(61) 



(62) 



(63) 



(64) 



This condition is not necessarily satisfied in the whole domain of convergence D. Therefore, together with the bound 
(|63p . this can be used to constrain modified gravity models. 



D. Retardation of light 

Let us study the gravitational time delay effect in light signals. For a light signal propagating from p to p, the 
integration of Eq. (f2"4")l with respect to p gives 



t = r s / dp 

J Pa 



Aq p 1 



-1/2 



(65) 



Expanding the integrand in e and assuming that both po and p are much greater than unity, the integral (|65[) is 
approximately given by 



t(Po,P) 



r s VPo - 1 



dp 



(p-l)Vp^pV 2 



p p [p 2 ( Cl +c 2 )-p 2 (2c 1+ c 2 )] + Cl pg 
(P 2 -P 2 f 2 



p+2 



(66) 

where p represents the position of the satellite. See Ref. [39| for the similar calculation about the gravitational 
time-delay effect. 
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The second term in Eq. (|66[) corresponds to the e contribution, t e . Under the condition p> po, the e contribution 
can be evaluated as 



P = l, 
P = 2, 
p = 3, 
P>4, 



r7 1 t e (p ,p) 



-- (ci + c 2 ) £( o 2 + i ep 2 [(c2 _ 3ci) + 2(ci _ Cs) m (2p/ Po )] , 

-g (ci +c 2 )ep 3 + - (ci - c 2 )£/3qP+ jce^, 

-i (ci + c 2 )ep A + \{ci- c 2 ) £ PoP 2 - ^ cepg [13 - 12 m(2p/p )] , 



1 Cl+C2 Fn p+i 

2 p+1 P 



1 Ci - C 2 2 p_j _ 1 

4 p-1 £PoP 4(p-3) 



ceplff 3 . 



where we have introduced the constant c defined as c\ when ci = c 2 . 

If p > 4 the time difference between two points p\ and p 2 coming from the e contribution is 

A , r s c l + c 2 / p+1 p +1n 



(67) 
(68) 
(69) 
(70) 

(71) 



The two contributions add because in the first integral one has propagation from the satellite to Sun (dtjdp < 0), and 
in the other one from Sun to Earth (dt/dp > 0). Since the standard GR contribution to Eq. (|6"6"|) is icR ~ r s ln[2p/po], 
the time difference between two points p\ and p 2 can be estimated as 



At 



GR 



r„ In 



4pip 2 



Pi 



The ratio among the two contributions is then given by 



At R 1 ci + c 2 p\ +1 + 



At GK 2 p+1 ln(4p 1( o 2 /p2) 



(72) 



(73) 



The bound regarding the ratio between the measured delay and the one predicted by GR comes from the Viking 
mission on Mars, which gives the result At exp /AiGR = 1-000 ± 0.001 [4(|. Hence this gives the bound 



At. 



At G R 



< 1(T 



Setting pi ss p2 ~ p, this condition translates into 

e|ci + c 2 | 



< 10 



-3 



2(p+l) In(2p/po) 

This is generally satisfied in the domain of convergence, as the logarithmic term grows slowly with p. 



(74) 



(75) 



E. Gravitational Redshift and Equivalence Principle 

Let us finally consider the gravitational redshift. In this case, for a light signal propagating at different heights r 
and r± , the ratio of the frequencies v and v± for corresponding heights is given by 



The e-dependent term is much smaller than the standard GR one in the D domain. Defining Av = v — v\, it then 
follows that 

Av e /v _ c x epp\(p v -p?) ^ 



Avqk/ v p- Pi 

The bound on [Av F /v)/ (Ah/ qr/^) comes from the experiment of an hydrogen-maser clock on a rocket launched to an 
altitude of about 10 7 m 4l|, which corresponds to Ai/ cxp /Ai/gr = 1 ± 0.0002. This leads to the following bound 

el <2xl0~ 4 . (78) 



Av GK /v 
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In the non-relativistic limit, the gravitational potential V can be identified as 500 = 1 + 2V, that is 



V 



1 



1 



GnMq 



1 

2 C1< 



2G N M G 



(79) 



This only depends on the mass M Q of Sun, not on the mass/properties of the test particle. Hence all test particles 
with same distance from the center will feel the same acceleration and the equivalence principle will not be violated. 

V. POWER-LAW f(G) MODEL 

The approach we have used so far works only if the iterative parameter e is much smaller than 1, and the method 
works better for smaller e. As a result this method cannot be evidently applied to all forms of R + f(Q). For example, 
let us consider the simple power-law case [42j 



f(G) = \VG* 



(80) 



such that f(G) is defined for all real values of k and Q . In fact, for the spherically symmetric spacetime, this Lagrangian 
will give rise to terms typically of order 



Qf,0- f = X(2k-l)eVG l 
Qs f,oo = 2\k(2k-l)e^/g. 



-, fc-i 



where 



-4fc 



— y/ G*/Gs 



l-Ak 



(81) 
(82) 

(83) 



The fact that e can be larger than one, implies that the corrections to the Schwarzschild metric may become large 
unless 1 — 4fc > 0, that is k < 1/4. The GR case corresponds to k = 1/2, in which the contribution of the GB term 
vanishes. If k 7^ 1/2 but close to it, then we still require that X\2k — l|/e < 1 in order to regard these terms as 
corrections to the equations of motion. In general, if k 7^ 1/2 and k > 1/4, the iterative method cannot be used to 
find approximate solutions to the equations of motion. In such cases the solutions need to be obtained by numerical 
integrations. 

If k < 1/4 then fgg is proportional to [{Gs/G) 2 ] 1 ~ k , so that this term blows up as G — > 0. In cosmological 
backgrounds the GB term G = 24H 2 (H 2 + H) changes sign from negative to positive during the transition from 
the matter era to the accelerated epoch [27| . This leads to the divergence of f t gg at Q = 0, which means that 
the pure power-law f(G) model is not cosmologically viable. We also note that the GB term inside and outside a 
spherically symmetric body (mass Mq and radius r Q ) with homogeneous density are given by G = —^(GnMq^/tq 
and G = 48(GjvAf Q ) 2 /r 6 , respectively. As we move from the interior to the exterior of the star the GB term also 
crosses from negative to positive. Although it is possible to derive iterative spherically symmetric solutions for 
k < 1/4 by using the expansion in terms of e, the power-law f(G) model is out of our interest because of the problems 
mentioned above. 



VI. MODEL A 

In this section we study the model (A) given in Eq. ([2]), i.e. 

f(G) = e [-a\^G~ s + \VG~ s g(x) 
where 

1,2 G 
g(x) = xarctanx ln(l +x ) and x = — — - 

2 Gs£ A 



(84) 



(85) 
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The Lagrangian is function of G/G*, but we choose to write it in this form so that the dependence on e becomes 
explicit. Since f gg > for positive A and G*, there is no singularity of this quantity unlike the power-law f(G) model. 
Let us discuss then the different contributions to the equations of motion (fTTj) . For the model (|84|) we have 



Gf.g -f = G- 



X 



<G s e 



arctanx + aXe^/Gs — XeyGsgix) 



i— 1 i— f G 2 
a Av^T+-A^:infl + ^ 



(86) 



which is of the order of eaX^fGl plus a logarithmic correction, 
motion, i.e. f^gg and f,ggg, can be written as follows 



The other terms that appear in the equations of 



f,gg 
f,ggg 



x 



fG~ s &g 

G 2 s 3 dx 2 



XVG~s d 3 g 
G 3 e 5 dx 3 



X\fG~ s £ 



G 2 e 4 



2XVGlGe 
(G 2 s 4 + G 2 ) 



(87) 
(88) 



which are both of the order of O(e) (as typically \G\ S> G s m the solar system). 
We write the metric in terms of the expansion parameter e: 



ds 2 = - 



1 - - + e0i (r) + e 2 fa{r) 
r 



di 2 



1 - — + efa(r)+s ^(r) 
r 



dr 2 +r 2 (d6 2 +sm 2 9dip 2 ). 



(89) 



Although the second-order contribution is not used in order to obtain experimental bounds on the model, we will 
evaluate it to check whether the terms in the series become smaller for higher orders of e and to verify that D 2 C D. 
Because of the Bianchi identities we can use only two equations, i.e. the 0-0 and the 1-1 equations, as the others are 
automatically satisfied. The second-order equations follow after solving the first-order equations. 



A. Spherically symmetric solutions and the domain of convergence 

Linearizing the 0-0 component of the modified Einstein equations at first-order in e, we obtain the differential 
equation for ipi in terms of p = r/r s : 

p-P + V'i = 32V3Xp 3 + 12V3Xp 2 ln(p) + (4 In e - 2a - 28)V3Xp 2 , (90) 
dp 

whose particular solution is 

V>i =8v / 3Ap 3 + 4\/3A i o 2 lnp+^\/3(21ne-a-16)A / 9 2 . (91) 

Here we have neglected the contribution coming from the homogeneous solution, as this would correspond to an order 
e renormalization contribution to the mass of the system. Although e 1 the term in lne only contributes by a 
factor of order 10 2 . Hence the first term on the r.h.s. of Eq. (f9"Tj) dominates over the last term. 
The 1-1 component of the Einstein equations gives the equation used to determine <fi as follows 

(P - P 2 )^ + 0i = 8V3Xp 4 - 2V3(10 + 6 lnp + 2 lne - a)Xp 3 - 2\/3(a - 6 Inp - 2 lne - 6)Xp 2 + ptpi . (92) 
dp 

Substituting the solution (|9"Tj) into Eq. (|92p . we get the following solution for fa: 

16 2 

fa = - — \/3A / 9 3 + -V3(4-a + 61n / 9 + 21ne)A/r 2 . (93) 
o o 

Since p > 1, the largest contributions to ipi and fa correspond to the ones proportional to p 3 , which are different 
from the Schwarzschild-de Sitter contribution (which grows as p 2 ). Hence the model ([2]) gives rise to the corrections 
larger than that in the cosmological constant case by a factor p. 
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In the Appendix we present the equations for the second-order quantities ip2 and 4*2- Discarding the homogeneous 
part of solutions of Eqs. (|A1[) and (|A2|) . we obtain 

V>2 = -512 A 2 p 7 - ^-A 2 (-3972 + 168 hie + 504 lnp-84a)p 6 - — A 2 (-322 lne + 161 a -966 lnp+ 1547) p 5 , (94) 

02 = A 2 p 7 - ^j-A 2 (-14a + 28 lne+ 1014 + 84 lnp)p 6 - — A 2 (294 hip - 399 - 49 a + 98 lne)p 5 , (95) 

where the dominant terms are the first terms in Eqs. (|94|) and (|95l) . 

So far we have found the solutions up to the order of e 2 . This expansion is meaningful if each term of the expansion is 
smaller than the previous one. For the solar system experiments, we will consider the case p> 1. This automatically 
imposes that the Schwarzschild contribution r s /r is smaller than the Minkowski value 1, which corresponds to the 
weak field approximation. The first-order correction in e is smaller than the Schwarzschild contribution if 

£|^i|<l/p, and e\<j)i\ < 1/p. (96) 

Since the dominant terms in ipi and </>i have the dependence Ap 3 , the conditions (|96p translate into 

(0 <|A£|- 1 / 4 . (97) 

Therefore, at first order, the domain of convergence is 

l«p«10 n |Ar 1/4 . (98) 

The consistency of this inequality requires that |A| <C 10 44 . The second-order terms in e can be neglected if the 
following conditions are satisfied in the same region 

e 2 -02 < eipi and e 2 </> 2 < efa . (99) 

Since the dominant contributions in -0 2 and 4>2 have the dependence A 2 p 7 , the conditions (|99[) arc equivalent to the 
requirement ([97]) . Hence our solutions using the expansion in e can be justified in the domain of convergence D. 



B. Solar system constraints 

We have seen that in the domain D of convergence the second-order term proportional to e 2 can be neglected. 
Moreover, in this domain, we can approximate the solutions further, keeping only the highest power in p = r/r s at 
the order e. Hence we have, for local gravity considerations, the metric (|12j) with 

i i a 

A=l — V3XEp 3 [l + 0{e,p- 1 )], (100) 

p 3 

B= 1- - + 8V3\sp 3 [l + 0(e 7 p- 1 )}. (101) 
P 

Therefore, in the domain of convergence, we will always regard the order e quantity as a correction to the Schwarzschild 
contribution, and for the order e we will always keep only the highest power in p. Since the solution given above is 
an approximate one valid in the domain of convergence, it does not make sense to see whether or not this metric is 
asymptotically flat: this would correspond to the solution in a region outside the domain of convergence. However, 
even having the real solution at hand (which is not the case), one should not trust it far away from Sun, as other 
forces would provide large contributions. 

In the following let us place constraints on the model parameter A by using a number of experimental bounds 
discussed in Sec. [IVJ 

• (A) Deflection of light 
The constraint (|37|) coming from the deflection of light gives 

2(p- l) £P 0PEarth < 10 ■ ( 102 ) 



For this model we have that 



16 

ci = -— \/3A, c 2 = 8V3A, p = 3. (103) 
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These numbers define the model and are the same for all the remaining constraints. The radius of Sun in units 
of the Schwarzschild radius is po = 2.35 x 10 5 . The distance of Earth from Sun in units of the Schwarzschild 
radius is pEarth = 5.08 x 10 7 . This translates into the following bound 

A < 1 x 10 15 , (104) 

where we have used the value (|9]) for e. 

• (B) Cassini experiment 

As we have already seen, the Cassini experiment places the bound y e < 10~ 14 , i.e. 

^Earth^^ ep Q (p P Ca L + PEarth) < 10"" ■ (105) 

The speed of Earth in units of the speed of light is t>Earth = 9.93 x 10 -5 , whereas the distance of Saturn in units 
of the Schwarzschild radius is pcass = 4.85 x 10 8 . This gives 

A < 2 x 10 12 . (106) 

• (C) Perihelion shift 

The bound (f6"3"|) coming from the shift of the perihelion of Earth leads to 

i |d| ep{p + 1)(1 + 2& art hK+ 2 th < \ ■ (107) 

For the eccentricity #Earth = 0.02 and the perihelion of Earth in units of the Schwarzschild radius pEarth = 
4.98 x 10 7 , we obtain the constraint 

A < 2 x 10 5 . (108) 
For Mercury, 5n g = 0.2 and pn g = 1-56 x 10 7 , the bound is slightly weaker, A < 5 x 10 5 . 

• (D) Retardation of light 

The bound (|74|) coming from the retardation of light, together with Eq. (f72|) . gives the following constraint 

1 g |Ci+g| p Mal , + p Eartl , < 1Q _3 ) (1Q9) 

2 p+1 ln(4p Ma rsPEarth/Po) 

which gives, with pMars = 7.71 x 10 7 in units of the Schwarzschild radius, 

A<5xl0 12 . (110) 

• (E) Gravitational redshift 

From the constraint |78|) coming from the gravitational redshift, together with Eq. (f77|) . it follows that 

|ci|£EarthP2Pl(P2 ~ Pi) <2x 1Q -4 ( m ) 



P2 - Pi 

Here £Earth ~ 10 -57 , pi is the radius of Earth in units of its Schwarzschild radius, i.e. p\ = 7.18 x 10 8 , and 
P2 = 1.84 x 10 9 is the distance of the experimental apparatus (for a height of 10 4 km). We then obtain the 
following bound 

A<3xl0 15 . (112) 



The tightest constraint on A comes from the perihelion shift experiment. This bound is weak so that the f(G) 
model ([2]) can be consistent with solar system constraints for a wide range of the model parameter. 
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VII. MODEL B 

Let us next proceed to the constraints on the model ([3]). This model can be written as 

f(G) = e [-a\y/g~ s + \y/g~.g(x)] , 

where 



(113) 



)(x) = x arctanx, 



and 



Q s e 2 



(114) 



We shall derive vacuum solutions for the spherically symmetric metric (|12[) by using the same expansion parameter e 
defined in ([8]). The term Qf.g — f in the equations of motion can be estimated as 



Gf,g - f = £^Qs 



£ 2 £ 4 + Q2 



Since Q z > G 2 e 4 , it follows that 



Of,g- f ~ e\^g s {a + \) + 0{e 5 



(115) 



(116) 



which works as a cosmological constant at lowest order. Note that we have f t gg ~ e 5 + 0(e 9 ) and f.ggg ~ s 5 + 0(e 9 ), 
so that these terms are higher than the linear order in e. These properties are different from those in the model ©. 

For the model ([3]) the dominant contribution to the Schwarzschild metric comes from the linear term in e and the 
next order corrections correspond to terms in e 5 . Hence we look for a metric of the form 



d^ 



1 



+ s<f>i(r) +s 5 ^> 2 (r) 



df 2 



1 



+ eipi(r) + s 5 tp 2 (r) dr 2 + r 2 (d0 2 + sin 2 6 d(p 2 ) 



(117) 



Linearizing the 00 component of the Einstein equation (|10p at first order in e, we obtain the following differential 
equation 



which has the particular solution 



p?i-^i + 2V3A(l + a)p 2 
dp 



A=-^X(l + a)p 2 



0, 



The 11 component of the linearized Einstein equation gives 



(P 2 - p)^ <h + ^A(c + l)p 3 - 2V3A(a + l)p 2 



whose particular solution is 



2V3 



A(l + a)p 2 



(118) 



(119) 



(120) 



(121) 



These expressions for the metric corrections represent the contribution of an effective cosmological constant (as in the 
case of Schwarzschild de Sitter metric). This solves the equations of motion up to the order e. 
Let us derive next-order solutions fair) and 02 {r). The differential equations for tp 2 and (j) 2 are 



n^P + V> 2 ~ 2V3A/9 14 (128p - 123) = . 
dp 



(p -pH — 

dp 

which have the following particular solutions 



lj'2 



2^3 



Ap 14 (80p 2 - 146p + 65) = 0. 



^Ap 14 (40 P -41) 
5 

2v / 3 , 14 
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Ap i4 (16 i o- 13) 



(122) 
(123) 

(124) 
(125) 
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In principle one can repeat this method order by order in e. It is then clear that the vacuum solution will not be 
strongly constrained by solar system experiments, as the first non-zero contribution will be that of a cosmological 
constant. To be more precise, the domain of convergence D is defined as 

eA(l + a)p 2 < 1/p < 1 , (126) 

which implies that 

£ -l/3 10 15 

l«p« 777: 77777 ~ 7777 777777. (127) 

[A(l + a)] 1 /3 [A(l + a)] 1 /3 

For the consistency of this inequality wc require that A(l + a) <C 10 45 . If A and a are of the order of unity, the domain 
D corresponds to the distance 10 5 cm < r < 10 20 cm. 

The contributions e 5 ip2 and s 5 (f>2 can be negligible relative to the first-order contributions provided that Ae 5 p 15 <C 
A(l + a)ep 2 , i.e. 

p< (l + a) 1/13 £- 4/13 . (128) 

If a = 0(1) the domain D 2 of convergence corresponds to r < 10 19 cm. One may wonder if this trend continues 
at next order, that is, e 6 . When the previous two contributions are coupled, this gives rise to terms of the order 
A 2 (l + a)e 6 . Introducing the corrections e e ip3 and e 6 03 to the metric and expanding the equations of motion at order 
e 6 , wc obtain the following differential equations 

P-P- + V J 3 - A 2 (a + 1) p 17 (608p - 1136) = , (129) 



dp 

A dp ™ " y " ' ^ ' 45 p 5 J 



, 9 A4>3 9, , 17 ( 9 1064 304\ 

(p 2 - p)^ - - \ 2 (a + I) p 17 [32p 2 + —p - — ) =0, (130) 



whose solutions are given by 



^-A 2 (l + a)p 17 ^P- 5 -f) , (131) 

3 = A 2 (l + a)p 17 (jP+ 1 ~§) ■ (132) 

Therefore the expansion is meaningful for A(l + a)e 6 p ls -^i s 5 p 15 , that is, D 3 = D C D 2 - This shows that the domain 
of convergence at the order e 6 coincides with 

Let us discuss solar system constraints on the model <j3j) by using the experimental bounds discussed in Sec. IIVI 
We will show that the strongest bound comes from the shift of the perihelion of Earth. 

• (A) Deflection of light 

For the model Q we have that 

ci = C2 = ~^-A(l + a) , p = 2. (133) 

Since c± = c 2l the constraint (|3T[> coming from the deflection of light is trivially satisfied. 

• (B) Cassini experiment 

Since c\ = c 2 at order e, the bound (|43|) of the Cassini experiment is fulfilled. 

• (C) Perihelion shift 

The bound coming from the shift of the perihelion of Earth corresponds to (|10T|) with ^Earth = 0.02 and 
PEarth = 4.98 x 10 7 . This leads to 

A(l + a) < 1 x 10 14 . (134) 
From the bound of Mercury we obtain a similar constraint. 
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• (D) Retardation of light 

Using the bound (|109p coming from the retardation of light with pMars = 7.71 x 10 7 , we obtain the following 
constraint 

A(l +a) < 6 x 10 20 . (135) 

• (E) Gravitational rcdshift 

Using the constraint ()111|) of the gravitational redshift with £Earth = 10 -57 , pi = 7.18 x 10 s , and pi = 1.84 x 10 9 , 
it follows that 

A(l +a) < 3 x 10 25 . (136) 
The bottom line is that the model B, having p = 2, is less constrained than the model A. 

VIII. CONCLUSIONS 

In this paper we have discussed solar system constraints on f(Q) gravity models that are cosmologically viable. These 
models give rise to power-law corrections of the form (r/r s ) p to the Schwarzschild metric, which are characterized by 
an expansion parameter e = \/G*/G s ~ 1CP 46 for Sun. The smallness of this parameter allows us to find approximate 
vacuum solutions in a spherically symmetric spacctimc. 

In order to confront f(G) models with a number of solar-system experiments, we have carried out general analysis 
for estimating their deviation from General Relativity These include the experiments such as deflection of light, 
Cassini tracking, perihelion shift of Earth, retardation of light, and gravitational redshift. The results we have derived 
can be generally applied to any modified gravity models which have power-law corrections to the Schwarzschild metric. 

The f(Q) models given in Eqs. ([2|) and ([3]) arc designed to give rise to a late-time cosmic acceleration preceded by a 
matter-dominated epoch. We find that these models can satisfy all of solar system constraints discussed in literature 
for a wide range of model parameters. For the model © there exists a logarithmic correction {X/2)^/Gl\n(G 2 /Q 2 ) to 
the Lagrangian density in the region of high density (Q 2 ^> Q 2 ). The tightest bound comes from the shift of perihelion 
of Earth, but the constraint on the parameter A is weak: A < 2 x 10 5 . In order to set stronger bounds on these 
theories, it is then necessary to have better measurement of the quadrupolc moment of Sun, as it affects the perihelion 
shift. For the model © the leading correction term to the Lagrangian density corresponds to (A / 'i)\/Ql Q 2 / 'Q 2 , whose 
effect is very small even compared to the model ([2]). Hence the model parameter for the model ([3]) is very weakly 
constrained: A(l + a) < 10 14 . 

The main reason why the f{Q) models can satisfy solar system constraints fairly easily is that, even in the vacuum 
spherically symmetric background, the Gauss-Bonnet scalar takes a non-vanishing value Q = 12r 2 /r 6 , where r s is the 
Schwarzschild radius. In the solar system the GB term is much larger than the cosmological value G* ~ Hi- Hence 
the inverse power-law terms such as (Q 2 /Q 2 ) n (n > 0) are strongly suppressed even for the vacuum solution. This 
property is different from f(R) gravity in which the Ricci scalar R vanishes in the vacuum spherically symmetric 
background. In this case the presence of non-relativistic matter with density p m leads to a non-vanishing Ricci 
scalar R approximately proportional to p m . Due to the existence of matter, the local gravity gravity constraints can 
be satisfied for viable f(R) models having the asymptotic behavior f(R) = R — XR* [l — (i? 2 /i? 2 )"] in the region 
R 2 3> R 2 - In f(G) gravity we have shown that the contribution of matter density p m to the GB term (~ p 2 ^) can 
be negligible relative to the vacuum contribution (~ r 2 /r 6 ) outside the area of Sun. Thus our analysis based on 
the vacuum spherically symmetric solution is reliable to discuss the compatibility of f(G) models with solar system 
experiments. 
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APPENDIX A: SECOND-ORDER EQUATIONS FOR THE MODEL A 

In this Appendix we present second-order equations in terms of the expansion parameter e for the model {2J . The 
second-order quantities tp 2 an d 4>2 m the metric (|89[) obey the following equations of motion 

(p 4 -p-3p 3 + 3 P 2 )^P- + (3p- 3p 2 - 1 + p 3 ) ip 2 + 4096 p 10 A 2 + (672 hip + 224 lne - 17488 - 112 a) X 2 p 9 
dp 

+ (-1040 lne + 520a - 3120 hip + 29472) X 2 p s + (-24448 + 1776 lne - 888 a + 5328 hip) X 2 p 7 

+ (664 a - 3984 In p + 9952 - 1328 lne) A 2 p 6 + (-1584 - 184 a + 368 lne + 1104 lnp) X 2 p 5 =0, (Al) 

and 

(-21 p 2 + 21 p 3 + 7 p - 7 p 4 ) -P- + (7 + 7 p 2 - 14 p) <j) 2 + 8512 p 10 A 2 

dp 

+ (-34976 + 1 12 a - 224 In e - 672 In p) X 2 p 9 + (168 In p - 28 a + 56 In e + 56260) A V 

+ (1344 In e + 4032 In p - 44440 - 672 a) A 2 p 7 + (- 1960 In e - 5880 In p + 980 a + 17 '444) A 2 p 6 

+ (2352 In p - 392 a + 784 In e - 2800) A 2 p 5 = . (A2) 

The solutions to this equation, discarding the homogeneous parts, are given by Eqs. ([94]) and (|95)) . 
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